Effective Langevin equations for the pair contact process with diffusion 



in 
o 
o 

(N 



(N 



43 

o 

CD 



CO 



I 

C 

O 

o 



> 

in 
o 
in 
o 



c 

o 
o 



X 



Ivan Dornic, 1 Hugues Chate, 1 and Miguel A. Munoz 2 

1 CEA - Service de Physique de I'Etat Condense, Centre d' Etudes de Saclay, 91191 Gif-sur-Yvette, France 
^Institute de Fisica Teorica y Computacional Carlos I, Universidad de Granada, 18071 Granada, Spain 

(Dated: February 2, 2008) 

We propose a system of coupled, real- valued, effective Langevin equations for the nonequilibrium 
phase transition exhibited by the pair contact process with diffusion (and similar triplet and quadru- 
plet, n-uplet, processes). A combination of analytical and numerical results demonstrate that these 
equations account for all known phenomenology in all physical dimensions, including estimates of 
critical exponents in agreement with those reported for the best-behaved microscopic models. We 
show in particular that the upper critical dimension of these n-uplet transitions is — , and — — 1 for 
their anisotropic (biased) versions. 
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Whereas there exists ample evidence of universality 
among out-of-equilibrium phase transitions, their full 
classification is by no means achieved. At equilib- 
rium, coarse-grained descriptions in terms of Ginzburg- 
Landau- Wilson free-energy functionals and their associ- 
ated Langevin equations encode in a systematic way the 
symmetries, conservation laws, and dimensionalities de- 
termining universality classes. Out of equilibrium, even 
the relevant ingredients are still debated, and contin- 
uous stochastic descriptions are often missing, so that 
the non-equilibrium counterpart of the celebrated taxon- 
omy by Hohenberg and Halperin remains a distant 
achievement. Here we are interested in irreversible phase 
transitions into fluctuation-less absorbing states. Such 
transitions abound in non-equilibrium physical phenom- 
ena such as epidemics, catalysis, synchronization, or self- 
organized criticality |2j • The prominent directed percola- 
tion (DP) class, which encompasses all phase transitions 
into a single absorbing state (without additional sym- 
metries, conservation laws, or disorder), well-established 
both theoretically and numerically, can be characterized 
by the following Langevin equation for a single coarse- 
grained, real density field p — p(x,t): 



ap — bp 2 + a p 2 r] 



(1) 



where rj = r)(x,t) is a Gaussian white noise whose influ- 
ence vanishes in the absorbing state (p = 0). In spite 
of some recent progress in the search for Langevin de- 
scriptions of other classes of absorbing phase transitions, 
the situation of the DP class remains exceptional. In- 
deed, even if Eq. can be derived rigorously for some 
reaction-diffusion models in this class, such Langevin 
equations sometimes simply do not exist. For example, 
for the simple annihilation process 2A — ^ it is well- 
known that one obtains "imaginary noise" (see below), 
which precludes a well-behaved Langevin equation 

Reaction-diffusion systems where creation and annihi- 
lation processes require pairs of particles (while isolated 
ones can diffuse but not react) are usually considered to 
represent the pair-contact process with diffusion (PCPD) 



class. Studies of various models in this (putative) class 
have produced a series of conflicting results and opinions 
(see |6j for a review and references therein). From the nu- 
merical side it is still debated whether in d = 1 the PCPD 
behaves as DP or not, as numerics are often plagued with 

long transients HUH E3, E3- what seems to be 
widely accepted is that the upper critical dimension d c is 
not 4 as in DP, but d c = 2 [H [l3| . It has also been 
shown numerically that when a bias (anisotropy) is in- 
troduced in the diffusion of isolated particles, the critical 
dimension is reduced to d c = 1 (at odds with what hap- 
pens for DP where biased diffusion plays no significant 
role) 0. At the analytical level Janssen et al. |l5j re- 
cently showed that a perturbative renormalization group 
analysis of the field theory derived from the microscopic 
reactions 2A — > 3 A, 2 A —*%/A (together with diffusion of 
isolated A particles) yields only unphysical fixed points 
and runaway trajectories. Generalizations of the PCPD 
requiring triplets or quadruplets of particles for the reac- 
tions (respectively called TCPD and QCPD for "triplet 
and quadruplet contact process with diffusion" ; n-CDP 
in general) have also been studied, with results at least 
as controversial as those for the PCPD 9, 16, 171. 



What makes these problems interesting beyond the 
usual skirmishes between specialists is that the ingredi- 
ent which seems to lead out of the DP class is the fact 
that 2 (n, in general) particles are needed for reactions 
to occur. That such a microscopic constraint, not related 
to any new symmetry nor conservation law, should deter- 
mine universal properties may indeed appear unaccept- 
able, hence the opinion, held by some authors 0, (13, EH , 
that the observed critical behavior will eventually turn 
up to be in the DP class. Elucidating the critical be- 
havior of the PCPD has thus become a milestone in the 
ongoing debate about universality out-of-equilibrium. 

In this Letter, we propose a system of coupled, real- 
valued, effective Langevin equations for the nonequilib- 
rium phase transition exhibited by the PCPD (and of the 
allied ti-CPD models). Analyses of these equations allow 
us: i) to show that d c — j- with a reduction by unity 



2 



when a biased diffusion is switched on, ii) to reproduce 
the correct (mean-field) critical exponents above d c in all 
cases, and iii) to estimate by direct numerical integration 
of these equations critical exponents in agreement with 
those of the best microscopic models in any dimension. 

A continuous description of the PCPD problem, if pos- 
sible at all, should be in terms of at least two fields as 
converging conclusions signal 0, 0|. Indeed, a two- 
field Langevin description has been successfully proposed 
for the related pair contact process without diffusion, the 
prototypical model with infinitely-many absorbing states 
[3j. As for that case, here a "pair-field" and a "singlet 
field" is the most natural choice. To build up a set of 
Langevin equations, we start from the following micro- 
scopic two-species reaction-diffusion model in the PCPD 
family Q: pair-activity is represented by usual DP-like 
reactions B — ► 2B, B ~ > while isolated particles diffuse 
and undergo pair annihilation 2 A — > 0. The two species 
are coupled cyclically: B — ► 2A (a "pair" separates into 
2 particles) and 2A — > B (2 particles recombine into a 
"pair"). Using well-trodden techniques (see, e.g., [lflj)- 
one arrives at the following coupled equations: 

d t ipB = D B V 2 ip B + a B ipB - b B ip B + c b^a + °Bi> B r\B 
d t ipA = D A V 2 4> A -b A ip A + c A tp B +io- A ip A r) A , (2) 

where all coefficients depend on microscopic reaction 
rates and diffusion constants, and r) A and r\ B are Gaus- 
sian white noises (including some cross correlations not 
specified here). The mapping through the two-species 
model disentangles the competition between the real, 
DP-like component of the noise, and the imaginary one 
(arising from the annihilation reaction 2 A — > 0) j2jj. But 
now both fields ipA and ip B are complex, and even though 
their noise-averages do correspond to the local densities 
Pa and p B of A and B particles, the same is not true for 
higher-order moments, impeding physical intuition. 

However, perusal of simulation results of PCPD-likc 
models suggests that annihilation is only predominant 
deep inside the absorbing phase, and not in the active one 
nor around the critical point. Indeed, slightly subcritical 
quenches reveal a quasi-exponential departure from the 
critical power-law scaling, with the typical annihilation 
power-law decay (tr 1 ! 2 for singlets and t~ 3 ^ 2 for pairs) 
setting up only at a much later time (see Fig. ^ and also 
0). Correspondingly, the variances of the noise terms 
— which can be estimated from measurements of block- 
spin like correlations of particle densities in microscopic 
models — are initially DP-like (i.e. positive and propor- 
tional to the particle densities) in the near-critical decay 
regime, then change sign (around t = 10 4 in Fig. ^ see 
also [lj|), and finally become negative (imaginary noise). 

On the basis of the above observations, and since our 
aim is to build an effective Langevin description of PCPD 
critical points, we simply discard the imaginary noise 
term (and its associated anticorrelations), and modify 
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FIG. 1: (Color online) Evolution, during a slightly subcritical 
quench, of the densities of isolated particles A and pair par- 
ticles B, together with the absolute value of the singlet-field 
noise variance (Spa) 2 measured as the connected correlation 
function of the local density of ^4s calculated between neigh- 
boring blocks of size 64. Two-species PCPD model described 
in text implemented as in 9] (d = 1, 2 23 sites, 16 runs . 

the nonlinear terms for the singlet-field so that it decays, 
for the space dimension of interest and when the pair- 
field is set to zero, with the appropriate decay at least 
on average |2l|. The abstract fields ipA,i>B in Eqs.J2J 
then become the real particle density fields pa and p B , 
governed by the following effective Langevin equations: 

d t p B = D B V 2 p B + ap B - bp B + cp a A + o B p\r] 
dtPA = D A \7 2 p A ~ epA + fPB, (3) 

which themselves change with dimension: for d = 1, one 
has a = 3 since the singlet-field must decay like \j\ft for 
pure annihilation, while for d > 2 the mean-field pure- 
annihilation decay l/t imposes a — 2. In the marginal 
case d = 2, the annihilation decay is lnt/i, and the in- 
criminated terms in Eqs.© have to be completed by 
logarithmic factors (±p A / In (1/ pa))- From these "bare" 
equations, one readily realizes that some extra terms — 
in particular a cross term ujpaP B in the p B equation or 
a noise term a 1 \fps~r\A in the pa equation — are gener- 
ated perturbatively once the effect of fluctuations is taken 
into account. Thus these terms should be incorporated 
from the beginning, and disregarding now all noises and 
derivatives, one arrives at the mean-field decay exponents 
#¥ F = — and 9^ F — 1, in agreement with the known re- 
suits for PCPD (where a — 2 in high dimensions) [6|- 

After having verified the validity of our Langevin equa- 
tions in high dimensions (mean-field) we now study their 
behavior in low dimensions by direct numerical integra- 
tion. For simplicity we first integrate Eqs.©, and af- 
terwards check that the extra perturbatively generated 
terms do not alter the critical behavior. Note that the 
equation for the pair-field p B is almost identical to the 
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DP one JIJ and, therefore, the accurate and efficient inte- 
gration scheme recently presented in Q can be employed. 
Varying the "temperature-like" parameter a, an absorb- 
ing phase transition is observed, with (apparently) clean 
scaling behavior at and around the critical point in all 
tested physical dimensions. As expected from our model- 
building, subcritical quenches show a quasi-exponential 
departure from the scaling regime, followed by a late 
annihilation- like decay (without local anti-correlations). 
In d = 1 critical quenches allow to estimate the decay ex- 
ponent 9 = (}/v\\ = 0.19(1) (Fig.©. Measuring steady- 
state activity above threshold determines j3 = 0.33(2), 
while finite-size lifetime statistics at criticality or seed 
simulations yield indifferently z = v_\_/v\\ = 1.64(3) (not 
shown). Even though the quality of these results, at par 
with those of the "best" microscopic models, leads to 
suggest asymptotic values different from that of the DP 
class in d = 1, we prefer in view of recent large-scale 
simulations in PCPD problems [l]), to avoid drawing a 
conclusion on whether the "true" asymptotic behavior is 
DP-likc or not, based solely on such a small quantita- 
tive difference (for the DP class 9 ~ 0.16, (3 ~ 0.28, and 
z ~ 1.58.) We have extensively checked that these results 
are robust against both changes in parameter values and 
inclusion of the extra perturbatively generated terms. 

In higher dimensions, the difference with DP transi- 
tions becomes qualitative: our simulations in d = 2 show 
that at threshold both fields decay with the same ex- 
ponent 9 d=2 ~ 0.50(2), the remaining curvature in log- 
log plots (not shown) indicating the presence of loga- 
rithmic behavior for the ratio (pb)/(pa) characteristic of 
the marginal dimension (in agreement with microscopic 
models). For d = 3, simulations of Eqs.© with a = 2 
show different scaling at threshold for the two fields: 
9 A = 3 ~ 0.50(5) and 9jf 3 ~ 1.0(1), values according with 
those found in microscopic models (Fig.©, and with the 
previously found mean-field predictions. This, and the 
behavior in d — 2 suggests that d c = 2. 

We now turn our attention to power-counting argu- 
ments aimed at identifying the critical dimension analyt- 
ically and determining the relevance or irrelevance of the 
different terms of our Langevin equations. As usual in 
problems with multiple fields, one has some freedom to 
perform naive scaling analyses. If we enforce, as usual, 
the coefficient of dtps to scale as a constant (fixing in 
this way the dominant time-scale) and the associated ac- 
tion to be dimensionless, one is ineluctably led to the 
conclusion that the noise term in the ps-equation is the 
dominant nonlinearity together with the saturation term 
~bp 2 B , and that both become relevant below d c = 4 as in 
DP. This being in contradiction with our findings, it is 
mandatory to scale in a different way. If we now impose 
the coefficient of the other time-derivative, dtPA, to scale 
as a constant, we conclude that the main nonlinearity is 
the (generated) noise term of the pa equation, propor- 
tional to y/pB, which becomes relevant below d c = A/ a. 
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FIG. 2: (Color online) Numerical integration of Eqs.Q in 
d = 1 {D A = 0, D B = 0.25, 6 = c = e = / = l, a% = 2, 
Ax — 1, At = 0.25, single run on system of size 2 23 ). De- 
cay of the (spatially-averaged) densities (pa) (top) and {pn) 
(bottom) during near-critical quenches from initial condition 
pA — 0, pB = 1 for a = 1.4442, 1.4443 and 1.4444. Inset: 
same multiplied by t 019 . 

Therefore d PCPB = 2 (as a = 2 in PCPD above its criti- 
cal dimension), in agreement with our numerics. One can 
thus consider that the (activity) pair- field ps is "slaved" 
to the singlet field pa (which is the leading one, control- 
ling the scaling). This viewpoint is further supported by 
the numerical observation that the decay-exponent 9b 
experiences a discontinuity at d c where it jumps from its 
mean-field value 9b — 1 above d c to 9b = 9a at or below. 
Following the above discussion, our final minimal set of 
effective equations yielding the correct scaling laws (as 
we have verified numerically) both when fluctuations are 
relevant and when they can be neglected reads: 

d t pB = D B V 2 p B + apB - bp 2 B + cp a A + ujp A p B 
d t p a = D A V 2 p A - ep a A + fp B + a A p\T}. (4) 

Note the change of perspective from Eqs.(©, which leads 
to consider the PCPD (and related models) as a sys- 
tem of annihilating random walks wandering between 
high-activity (pair) clusters, themselves of course self- 
consistcntly determined by the complex interplay of the 
two modes of the dynamics. These clusters act not only 
as non-trivial boundaries for the isolated particles, but 
also as fluctuating sources. Since fluctuation effects for 
2A — * (or 2A — > A) in the presence of even a constant 
source term are known to be relevant up to and includ- 
ing d = 2 \22\ . this provides, in our opinion, an intuitive 
mechanism which can explain, among other things, why 
d c = 2 for the PCPD and also why the dynamical expo- 
nent is not that of simple random walks |23| . 

The two-species/two-fields description put forward 
above is easily extended to the TCPD and QCPD cases 
EJ fivt - by replacing the pair field by a triplet or 
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(a) (b) 




FIG. 3: (Color online) Decay of singlet and pair densities 
during near-critical quenches for d = 3 PCPD models, (a) 
Reaction-diffusion model 2A — > 3^4, 2A — > implemented as 
in 9J; (b) Integration of the Langevin Eqs. © with a = 2. 



quadruplet field: the interpretation of the fields changes 
but the governing equations |j2J remain identical, with an 
exponent a chosen according to the corresponding be- 
havior of the pure annihilation processes (3A — > or 
— ► 0). For 3A — > 0, the particle density decays as 
•y/ln i/t in rf = 1, which can be accounted for by a = 3 
with a logarithmic correction, while for d > 1 mean-field 
decay (corresponding to a — 3) emerges. A strik- 

ing consequence is that the TCPD class should not exist 
per se in d = 1, since its effective Langevin equations 
are identical to those of the PCPD class (the logarithmic 
term plays no role on asymptotic scaling). Moreover, 
following the power-counting analysis sketched above, 
^tcpd _ 4 a g reemen t w ith with the mean- 

field exponents 9 a = \ and 9 b = 1. Now for quadruplet 
models, the pure- annihilation mean-field decay holds 
in all physical dimensions and thus a = 4 always, lead- 
ing to d^ CPD = 1. Indeed, in d = 1 both fields decay 
with exponent 9 = i with strong log corrections (most 
likely ( Pa )/(pb) ~ (mt) 1/3 ), while for d > 1, 9 A = \ 
and 9b — 1 (mean-field results). All these results (in 
disagreement with some previously published ones) have 
been verified by simulations of the Langevin equations 
and of microscopic two-species reaction-diffusion models, 
and will be reported elsewhere. 

Finally, let us discuss the role of biased diffusion 
(anisotropy) for the isolated particles within our picture 
|l4j . Starting either directly from microscopic models, 
or from symmetry considerations, one realizes that an 
additional Vp^-term needs to be introduced in the p A - 
equation. A power counting analysis, analogous to the 
one discussed above, then reveals that d c is reduced by 
one unit, so d c = 1 for a = 2 (biased-PCPD) and it is 
below d = 1 for larger values of a: the biased TCPD and 
QCPD are expected to exhibit mean-field behavior in any 
physical dimension. All these results are in agreement 
with known ones [l4T | and have been verified numerically. 

Summing up, we have proposed effective coupled 
Langevin equations governing the densities of a n-uplet 



field and a singlet field in reaction-diffusion problems 
where at least n particles are required for creation or 
annihilation to occur. A combination of numerical and 
analytical arguments have shown that they do repro- 
duce the behavior of microscopic models, and thus pro- 
vide a sound framework to further study these problems. 
Such future work should concentrate on a renormaliza- 
tion group treatment of these equations, possibly in a 
non-perturbative approach. Due to the closeness of esti- 
mated exponent values for the PCPD case in d — 1 with 
those of the DP class, this appears as the only solution 
to put an end to the crucial question of whether such 
problems, below their upper critical dimension, exhibit 
genuinely novel critical behavior or not. 

M. A. M. acknowledges financial support from the 
Spanish MCyT (FEDER) under project BFM2001-2841. 
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